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We revisit the scaling relationships between the dark matter mass and observed X- 
ray luminosity and temperature of galaxy clusters and groups in the local Universe. 
Specifically, we compare recent observations with analytic models of the intracluster 
medium in which the gas entropy distribution has been shifted by a variable amount, 
K , to investigate the origin of the scatter in these scaling relations, and its influ- 
ence on the luminosity and temperature functions. We find that variations in halo 
concentration or formation epoch (which might determine the time available for low 
' entropy gas to cool out) are insufficient to explain the amount of scatter in the mass- 

^ \ luminosity relation. Instead, a range of entropy floors at a fixed halo mass, spanning 

& ■ approximately ~ 50 keV cm 2 to ~ 700 keV cm 2 , is required to match the data. This 

range is likely related to the variance in heating and/or cooling efficiency from halo to 
halo. We demonstrate that these models are consistent with the observed temperature 
■ and luminosity functions of clusters, with a normalization of erg ~ 0.8 in agreement 

5— i ' with WMAP measurements (for h = 0.7 and fl m = 0.3); in particular the scatter in 

the mass-luminosity relation has an important influence on the shape of the lumi- 
nosity function, and must be accounted for to provide a consistent result. Finally, we 
present predictions for the redshift evolution of these scaling relations and luminos- 
ity/temperature functions. Comparison with recent data at z < 0.7 shows reasonable 
agreement with a model that assumes a median entropy floor of K a — 200 keV cm 2 . 
When observations are extended to group scales (kT < IkeV) , this evolution will have 
the potential to discriminate between an entropy floor that is independent of redshift 
(for example, in a preheating scenario) and one that depends on the cooling time of 
the halo. 
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1 INTRODUCTION 

The X-ray properties of clusters are tracers of both the 
gravitational potential and the thermodynamic history of 
the gas. Since the mass is likely dominated by collision- 
less dark matter, for wh ich we have a well-developed theory 
fe.g. lEvrard et al.ll2002n . we might hope to learn about the 
relevant baryonic physics through detailed X-ray observa- 
tions. In particular, clusters obey fairly well defined scaling 
relations between mass and X-ray temperature (M-T) and 
mass and X-ray luminosity (M-L). Although the slope and 



normalization of these relations rel ative to model predic- 
tions have been studied extensively llEdge fc Stewartl| l99l|: 
MarkevitcfJl998l:lHorner et alll99alNevaIainen et all200C : 
Arnaud et alJ l2005'). litt le attention has been paid to their 
intrinsic scatter (b ut see lRowlev et al.l2004l : lKav et al.l2004l : 
ISmith et alj Eiol. 

The presence of scatter in the M-L and M-T relations 
suggests intrinsic variations in the structure of clusters, 
which may be due either to variations in the dark matter 
distribution of the halos themselves, and/or in the gas prop- 
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erties. In particular, the scatter is likely to be driven by vari- 
ations in the core properties of clusters, since the scatter is 
known to be significantly reduc ed if the central r egions are 
excluded from the analysis fe.g. iMarkevitcbl ll 19981 , 

In the case of the dark matter component of relaxed 
clusters, there are indications that either (or both) the halo 
concentration parameter and the inner slope of the halo pro- 
file varies from cluster to cluster, perhaps due to cosmic 
environmental effects such as the extent of tidal torquing 
the dark matter experienc es during c o llapse or the merger 
history of the halo fe.g. Ijind l20o3 iBullock et all EoOll 
IWechsler et, al.l2002l : lzhao et alJ200alWilliams et, alJ20f)4l . 
The amount of substructure and dynamical state of the dark 
matter will also vary and depend upon the merging history 
of the halo. 

There are also reasons to expect substantial variation 
in the gas properties of clusters, independently of their dark 
matter distribution. Models that seek to account for the 
mean M-L and M-T relations require some form of en- 
tropy modification i n the central reg ions of the systems 
jThomas et al.ll2002l : fviana et al.l l2003'l. Pure heating mod- 
els have been very successful in explaining the average trends 
of these scaling relations, especially when the heating tar- 
gets the lowest en tropy gas jKaiseJl99lHBalogh et al.ll999l 
iBabul et aDl200l hereafter BBLP). It seems likely that the 
efficiency of whatever physical mechanism is responsible for 
the heating (e.g. heat transport, AGN energy injection, etc) 
will vary from cluster to cluster. In this case, one expects 
both cooling and the feedback it triggers to eliminate gas be- 
low a n entropy threshold that depends on halo mass and red- 
shift llVoit fc Brvarj200ltlVoit et al.l20ol hereafter VBBB), 
although some gas may exist below this threshold if it is in 
the process of cooling out. In a realistic cooling model (e.g. 
iMcCarthv et al.l2004h . scatter can be introduced by appeal- 
ing to a range in time available for cooling in each halo. 
Despite this, McCarthy et al. show that such a range alone 
cannot account for the scatter in the M-L relation; it is also 
necessary to introduce variations in an initial heating level. 
However, they do not consider the effect of variations associ- 
ated with the underlying dark matter potential, which will 
contribute some scatter independently of the gas entropy 
distribution. 

Shorter-lived changes to the equilibrium temperature 
and luminosity of a cluster may also be associated with 
merger events (e.g. IRitc hie fc Thomasll2002l : iRandall et alJ 
120021: iRowlev et al-lhoolT Departures from equilibrium in 
the potential can change the luminosity, but ha ve little ef- 
fect on the gas temperature iRowlev et alj l2004'l. However, 
the shocks associated with mergers can have a temporary 
but significant influ ence on both the tempe r ature and lu- 
mino sity of the gas (IRitchie fc Thomas! I2OO2I IRowlev et~aH 
12004) . 

Recently, semi- and fully-numerical simulations which 
include both cooling and feedback from star formation have 
been shown to produce clusters with X-ray properties that 
scale with mass in a way that is in reasonable agree- 
ment with the observations fe.g. [Mjianwonj^^dJ |20M ; 
Thomas et al.ll2002l: IrBorgani et al]|2002l IViana et all 2003 : 
Kav et al] l2004t IRowlev et alJ 120041 : iBorgani et al.l I2OOE : 
Qstriker et aljl2005l) . [n this paper we will re-examine the 
observed M-L and M-T relations, focusing on the scatter 
in these relations and how it relates to the expected varia- 



tion in the underlying physical processes. The data we use 
are uncorrected for any cooling-core component, since it is 
precisely this core region that interests us. We will compare 
these observations with analytic, hydrostatic models which 
allow modifications to the entropy distribution of the gas 
(BBLP, VBBB), to determine the range of model parame- 
ters that are required to reproduce the observed scatter. Al- 
though the hydrostatic nature of these models means they 
are not ideally suited to explore in detail the effects of active 
cooling or heating in clusters, the range of model parame- 
ters required to match the data can be related indirectly to 
these processes. We will also make a self-consistent compar- 
ison with the temperature and luminosity functions, which 
provide an independent test of t he models ass uming the dark 
matter mass function is known fevrard et all2002T) . Finally, 
we will present the redshift evolution of all these observable 
quantities to put further constraints on the model parame- 
ters. 

Throughout this paper we use a cosmology with f2 m = 
0.3, tt A = 0.7 and H a = 70 km s" 1 Mpc" 1 . The theoretical 
models are described in § D2J and the comparison with ob- 
served X-ray properties is presented in § |H] Predictions for 
the evolution of these models, and some comparison with 
early data, are given in §|S] We summarize our conclusions 
and discuss the implications and limitations of our findings, 

M0 



2 THE MODELS 

2.1 Dark matter profile shapes 

The average theoretical shapes of dark matter halos are well 
motivated by N-body simulations to have a form given by 

p oc r~ ni (1 + c 2 oor)~ n2 (!) 

iNavarro et al.l ll996. NFW), where m, n%, and C200 are fit- 
ting parameters. The primary determinant of a halo's struc- 
ture is its virial mass, M v j r . This is commonly defined, us- 
ing spherical collapse models, as the mass within a fixed 
overdensity A that depends on cosmology a nd redshift; for 
AC DM, A - 100 at z = <Eke et al.lll996f) . However, ob- 
servations are more typically made at A = 200, 500 or larger. 
The radius and mass corresponding to the overdensity A will 
be denoted Ra and Ma, respectively. 

High resolution simulations show that relaxed clusters 
have remarkably uniform profiles. We will assume an NFW 
profile with ni = 1 and ri2 = 2 as our fiducial model. The 
range of v alues of these parameters reported i n the liter- 
ature fe.g. iMoore et alJll99Sl: iLewis et alll2000f) appear to 
be mostly due to differences in resolution and the fact that 
the fitting formula is not a perfect description of the profile 
(Hav ashi et alj |2004). Within a given simulation, an indica- 
tion of the amount of variation in halo shapes is best given 
by the distribution of concentration parameters, 0200- This 
parameter has a systematic dependence on mass and red- 
shift fe.g.lEke et al.|ll99l IBullock et alll200ll : |Powerll2003l : 
IWechsler e^L|[ 2002T) : we will take the parameterization of 
lEke et al] B2001J) . assuming as = 0.8. Most importantly for 
our purposes, C200 shows conside rable scatter at fixed mass 
and redshift; iDolag et al.l (l2004) have recently shown that 
the distribution of C200 values in simulated clusters (selected 
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Figure 1. Top panel: The relation between C200 (concentra- 
tion) and M200 as a function of Af2oo is shown as the shaded 
region, where the distribution is due to the 1 — a (heavy shad- 
ing) and 3 — a (lig ht shading) distribution of concentrations from 
Dolag et al. (2004). The horizontal line shows a fixed concentra- 
tion of 4, for reference. Bottom panel: The relation between 
M500 and A/200 as a function of M200 is shown as the shaded 
region, assuming the same distribution of concentrations above. 
The horizontal, solid line shows the result for a model with a fixed 
concentration, C200 = 4. 

only based on an overdensity criterion) is approximately log- 
normal with a width that is nearly independent of mass. 
This scatter is at least partly due to the presence of sub- 
structure, triaxiality, and departure from equilibrium in the 
sample of simulated clusters chosen from the simulations. 
We will therefore consider a range of concentrations corre- 
sponding to the ±3<T range predicted from this distribution. 
Although different values of as will change the value of C200 
at fixed mass by a small amount, it will not have an impor- 
tant effect our discussion in this paper, which is based on 
the variation in C200 and not its absolute value. 

In Figure we show how the concentration and its scat- 
ter depend on dark halo mass M200 in our model. The aver- 
age concentration and its scatter both decline with increas- 
ing mass. In the bottom panel we show how this affects the 
ratio A^oo/A/soo- The scatter in concentration at a fixed 
mass corresponds to a ~ 10 per cent scatter in A^oo/AZsoo- 
The dependence of concentration on mass means the slope of 
measured correlations between observables (like X-ray tem- 
perature and luminosity) and mass will depend on which 
mass is used. Throughout this paper we will present our re- 
sults as a function of A/200 , and this relation can be used to 
deduce the corresponding value of A/500 • 



2.2 Entropy distributions 

The shape of the gas profile in a cluster of given mass is 
determined by the entropy distribution of that gas, which 
is sensitive to its thermodynamic history. We adopt the 



common (e.g. iPonman et allll999l) redefinition of entropy 
as K = kT e rie 2 ^ 3 , where T e and n e are the electron tem- 
perature and density, respectively. This is related to the 
thermodynamic entropy by a logarithm and an additive 
constant, and is given in units of keV cm 2 . We will usu- 
ally quote this quantity in dimensionless units, relative to 
Kiqo = lOOkeV cm 2 . 

We use the formalism of VBBB to compute the hy- 
drostatic equilibrium gas distributions under different as- 
sumptions about the dark matter potential and the ther- 
modynamic history of the gas. Once the halo potential is 
specified, the entropy distribution and appropriate bound- 
ary conditions are all that are required to fully describe the 
gas density and temperature profiles. We start with an initial 
profile where the gas density traces the dark matter density, 
and solve for the temperature profile needed to satisfy the 
assumption of hydrostatic equilibrium. This is a good ap- 
proximation at large radii, where the entropy scales approx- 
imately as K oc r 1 ' , as found from analytic modelling and 
numerical simulations jLewis et al.ll2000t Irozzi fc Normanl 
l200ltlVoit et al.ll2003l: iMcCarthv et alJl2005l) . The gas dis- 
tribution at smaller radii is dominated by the entropy mod- 
ifications that we discuss below, so the choice of initial pro- 
file is relatively unimportant. The normalization of this ini- 
tial, unmodified model, is chosen by assuming the gas frac- 
tion within -R200 is equal to the global baryon fraction of 
12.9 per cent (for = 0.7 and ft = 0.3. iBurles et alJl200ll : 
ISpergel et aT]l2003D . 

We will explore modifications to this default profile, in 
the form of the shifted-entropy models of VBBB. In this 
case, the entropy distribution is shifted by an additive con- 
stant, which provides a good approximation to pre-heated 
models. This entropy shifting causes the gas to expand be- 
yond -R200; following VBBB we therefore choose as our pres- 
sure boundary condition the accretion pressure at the max- 
imum radial extent of the gas. The amount by which the 
entropy distribution is shifted will be left as a free parame- 
ter; for our base model we use K a = 2/\ioo, which is known 
to provide a reasona ble match to the medi an global scaling 
relations of clusters iMcCarthv et al.l l2004). 

We will also consider the case where K is set by the 
entropy of gas that can efficiently cool in time t. This en- 
tropy, which we will call K coo \, depends on mass and redshift 
(Voit & Brvan 2001), but also on the time available for cool- 
ing (IMcCarthv et alJ 120041) . Figure [5] shows the maximum 
value of K predicted in this cooling-based model, calcu- 
lated assuming the gas can cool for a Hubble time. The 
mean value of K coo \ decreases with increasing redshift due 
primarily to the decrease in time available for cooling. Of 
course, simply shifting the entropy distribution by a con- 
stant value is not an accurate representation of the effects 
of radiative cooling. In reality, an inner entropy gradient 
(rather than a floor) may be e xpected, with some gas below 
the entropy threshold (VBBB; IMcCarthv et alll2004l) . This 
can be particularly important for clusters, where the central 
entropy in a more realistic model can drop well below the 
value of if cool- Furthermore, the time available for cooling 
will likely vary from cluster to cluster, but must always be 
less than a Hubble time. Both of these effects tend to reduce 
the value of K a \ thus our prediction in Figure [5] represents 
a strict maximum value for this quantity. 

It is important to note that the entropy modifications 
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Figure 2. The entropy K coo \ corresponding to the maximum 
entropy for gas that can cool in a Hubble time, as a function of 
halo virial temperature and redshift. The horizontal, dashed line 
shows the K coo [ = 2A"ioo line for reference. 



that we explore in the present study are all confined to the 
cluster core. Therefore, our approach is to effectively bracket 
the range of central entropy levels required for the heating- 
based and cooling-based models to explain the scatter in 
the observational data. Both models implicitly assume that 
the entropy distribution at large radii is essentially identi- 
cal to that found in clusters formed in cosmological numer- 
ical simulations. Gas at large distance from the centre is 
not easily affected by cooling or non-gravita tional heating 
processes once the cluster is assembled (e.g. iBorgani et alJ 
120051: lOstriker et al-lEoO?! . However, heating of the gas be- 
fore it is accreted into a cluster can smooth out the density 
distribution of infalling gas, and this increases the entropy 
jump at the accretion shock llVoit et all2003HBorgani et alJ 
2005). In this case, the entropy distribution at large radii 
will be larger than that we have assumed, and a lower value 
of central entropy will be required to explain the observa- 
tions. This underscores the need for detailed comparisons 
of theoretical models to the spatially-resolved entropy dis- 
tributions of large, representative samples of clusters. At 
present, however, only a relatively small number of clus- 
ters have accurately-determined entropy profiles from Chan- 
dra and XMM-Newton data. Analysis of this small dataset 
seems to confirm that for high mass clusters (kT <; 4 
keV) the entropy distribution at large radii does indeed 
trace the spatial distribution and normalizatio n predicted 
by 'a diabatic' hydrodynamic simulations (iMcCarthv et al] 
2004). This conclusion is strengthened by the fact that for 
high mass clusters there also appears to be excellent agree- 
ment between the observed projected temperature profiles 
at large rad ii and those predicted by hydrodynamic simu- 
lations (e.g..lDe Grandi fe MolendJl200a lLoken et al.ll2002t 
IVikhlinin et alJl200rjl . Thus, our estimates of K a for such 
systems should be robust. However, for cooler systems, there 
are preliminary indications of excess entropy at large radii 



fe.g.. |Ponman et al-lEooS iPratt fc Arnaudll2005l) . As such, 
our estimates of K for low temperature systems should be 
treated with some caution. In this paper, our conclusions 
rest primarily on the data for the high mass clusters. 

2.3 Prediction of observable quantities 

To compute X-ray observables fr om the analytic gas pro files 
we use the cooling functions of iRavmond et al for 
gas with one third solar metallicity. To avoid the need to 
make bolometric corrections to the data, which depend on 
an accurate measurement of the gas temperature, we com- 
pute the model luminosities within the observed ROSAT 
energy bands 0.1-2.4 keV and 0.5-2.0 keV. Total luminosi- 
ties are obtained by integrating out to a minimum surface 
brightness of 1 x 10 -15 ergs s" 1 cm" 2 arcmin -2 , similar 
to that of the WARPS survey dScharf et al.lll997f) . unless 
stated otherwise 1 . The sensitivity of our results to this limit 
are explicitly noted. 

The model temperatures we compute are emission- 
weighted by the 0.1-2.4 keV luminosity, again excluding re- 
gions below the WARPS surface brightness limit; however 
the choice of energy band and surface brightness limit have 
a negligible effect on the calculated temperatures for our 
purposes. Simulations suggest that spectral temperatures, 
as measured observationally, can be 1 0-20 per cent h igher 
than emission-weighted temperatures |Mathiesen fc Evrardl 
120011 iRasia et alJl2005t IVikhlininll2005l) . For relaxed clus- 
ters, the difference is probably closer to the low er end 
(~ 10 per cent) of this range (Rasi a~et alJ 1200*51) . How- 
ever, this systematic error is not of major concern for the 
present study, as the statistical errors associated with the 
observed temperature, mass, and luminosity of a given clus- 
ter are typically twic e as large as this. For example, in 
IMcCarthv et al.N2004h we derived similar constraints on the 
parameters of heating and cooling models of the intraclus- 
ter medium (ICM) from independent analyses of the M-L 
and luminosity-temperature (L-T) relations, which suggests 
that small systematic errors in temperature measurements 
do not have a noticeable influence on our results. 

Our pr imary source of data is th e HIFLUGCS clus- 
ter sample llReiprich fc Bohringerll2002l) . This survey is an 
X-ray flux-limited sample of nearby clusters based on the 
ROSAT All Sky Survey, in which cluster masses are deter- 
mined from the density profiles, assuming hydrostatic equi- 
librium and isothermal temperature profiles. It is common 
practise when considering X-ray scaling relations to use tem- 
peratures and luminosities that are corrected for a cooling 
flow component, and it is known that this r educes the scat- 
ter in these relations fe.g. lMarkevitchll99ctl . However, since 
it is precisely this scatter that is the focus of our study, we 
will present all the data as observed, without this correc- 
tion. We take these raw cluster luminosities, in the 0.1-2.4 

1 We note that Scha rf et alJ il997f) attempt to correct for flux 
below the surface brightness limit. However, the flux correction 
to the data (typically a factor ~ 1.4) is an underestimate at the 
lowest luminosities, since it assumes a surface brightness profile 
slope of /3=0.67, appropriate for high and moderate luminosity 
clusters but not for low luminosity groups. Thus, the lowest lu- 
minosity data points may still systematically underestimate the 
luminosities by a factor ^2. 
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Figure 3. The relation between X-ray luminosity in the 0.1-2.4 
keV band and halo mass M2oo- The filled circles are local (z < 
0.2) data from the HIFLUGCS sample iRcip rich fc Bohringerl 
2002). 1— ct error bars on the masses are shown; only clusters 
with relative errors of < 50 per cent are plotted. The shaded 
region shows the shifted entropy model with Ko = 2i\ioo and 
a lcr (heavy shaded region) and 3a (lighter shaded) range of 
halo concentrations. The dashed line represents the model with 
Ko = K coo i, using the most probable value of C200 at each mass. 

keV band, from the lReiprich fc Bohringerl l|2002l) catalogue. 
Since many of the temperatures and masses in this cata- 
logue have been subjected to a significant cooling-flow cor- 
rection, we will only keep the ~ 80 clusters (out of 106) 
for which u ncorrected tem peratures are available, from the 
catalogue of lHornerl i200lfl . The dynamical mass estimates, 
A/200 and M500, are derived from the gas temperature, as- 
suming hydrostatic equilibrium. To be fully consistent, we 
make a small adjustment to these masses (M oc T) so they 
agree with the original (uncorrected) temperatures. 



3 SCATTER IN THE M-T AND M-L SCALING 
RELATIONS 

3.1 The M-L relation 

In Fig ure E3 we present the M-L relation from the data of 
iReiprich fc Bohringerl (l2002fl . excluding those clusters with 
mass uncertainties greater than 50 per cent. Note that the 
mass we plot here is JV/2oo- Observationally, M500 can be 
more precisely determined, and this c an reduce the purel y 
observational scatter in this Figure (iRowlev et al.ll2004) . 
However, when we construct the temperature and luminos- 
ity function in § 0] wc will have to use M200 , since this was 
the mass us ed to derive the m ass function from numerical 
simulations jEvrard et al.ll2002T) . For consistency, therefore, 
we have used M200 throughout the paper. We show both 
observed and model luminosities in the 0.1-2.4 keV band, 
to minimize errors in bolometric corrections to the obser- 
vations. The limiting surface brightness cut applied to the 



models (see § results in a significant reduction of lumi- 
nosity only below M200 ~ 1O 13 ' 5 M0. The shifted-entropy 
model reproduces the slope and normalization of the relation 
well. However, the scatter in the data is much larger than 
expected from th e observational uncertainties, as has been 
noted before (e.g. lFabian et al]|l994l ; lMarkevitchlll998l) . 

Our calculations show that the observed scatter is much 
larger than can be expected from variations in halo concen- 
tration. This is shown by the shaded region in Figure |3] 
which represents the 1— and 3 — a range of luminosities 
predicted at a given mass, from the dispersion in simulated 
cluster concentrations alone. Approximately 25 per cent of 
the observed clusters lie well outside the 3<r range resulting 
from variations in halo concentration. Even though our mod- 
els assume spherical, smooth, virialized halos, these effects 
are partly accounted for by the variation in concentration 
parameter, which is determined from simulated clusters that 
are clumpy, non-spherical and in a variety dynamical states. 
Furthermore, the predicted scatter in the scaling relations 
is not significantly larger in models which consider more re- 
alisti c potential shapes iRowlev et al 1 l2004l : lOstriker et alJ 
2005). It is therefore unlikely that the observed scatter in 
this relation can be entirely attributed to variations in the 
shape of the dark matter distributions. 

Some of the scatter in the M-L relation could be due 
to short-timescale events like mergers, which cause changes 
in both the luminosity and temperature of the gas (e.g. 
[ Ritchie fc Thomas! 12002 : iRandall et ai1l2002l ; IRowlev et ail 
l20^) ~lkowever. it has been shown that these changes tend 
to move galaxies along the L-T relation, and do not con- 
tribute significantly to its scatter; since the observational 
scatter in the L-T relation is compara ble in magnitude to 
that in the M-L relation of Figure H ( F abian et al 1 ll994l ; 
lMarkevitcf]|l998l) . there must be an important source of 
scatter other than mergers. Furthermore, as major merg- 
ers in massive clusters are expe cted to have been relatively 
rare in the past ~ 2 Gyr (e.g. iKauffmann fc White! Il993l : 
lLacev fc Colelll994l) . and the luminosity and temperature 
boosts typically last for ~ 0.5 Gyr or less following a ma- 
jor merger, these events are relatively rare and unlikely to 
be responsible for all the observed scatter in an unbiased 
cluster sample. 

On the other hand, the X-ray luminosity of a cluster 
of given mass is very sensitive to the entropy floor level, 
as shown in Figure 01 A range of entropy levels 0.5-5 -ft'ioo 
approximately covers the scatter in the observations, at high 
masses. For a few clusters even higher levels of K ~ 7A"ioo 
are required to match their low luminosities, something that 
was also observed in our earlier comparison with Sunyaev- 
Zeldovich measurements jMcCarthv et alJl2003l) . This is a 
remarkably large range in central entropy levels, and thus 
substantial variations in heating or cooling efficiency must 
exist from cluster to cluster, even if some of the observed 
scatter can be attributed to substructure, departures from 
equilibrium, or merger-induced shocks. 

In the cooling-based model, the entropy floor is related 
to the cooling time (as in Figure |2J . This model predicts 
that Ko depends on mass, increasing by a factor of ~ 5 from 
~ O.8-K100 to ~ 4Kioo over the temperature range of inter- 
est. The resulting M-L relation is shown as the dashed line 
in Figure^] the slope here is steeper than for any of the fixed 
entropy-floor models shown in Figure^] Note, however, that 
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Figure 4. The mass-luminosity relation, with data as in Fig- 
ure|3] and five models of fixed concentration but different entropy 
thresholds, as labelled. 



while K may vary by a factor of ~ 5 over this mass range, a 
similar or larger range of K is required to explain the distri- 
bution of luminosities at fixed mass. Although the cooling 
model can easily accommodate scatter toward more lumi- 
nous clusters (since -ftT C ooi is only the maximum entropy of 
gas that can cool), there is no simple mechanism to account 
for the scatter of clusters toward lower luminosities. There- 
fore the simple interpretation that the entropy threshold is 
due solely to the cooling of low entro py gas is not likely to be 
corr ect, a conclusion also reached bv lMcCarthv et "all J2004T) 
and lBorgani et ail j2005tl using more sophisticated models 
that account for the presence of gas that cools below the 
K coo i threshold. 



3.2 The M-T relation 

In Figure we compare the m ass-temperature relation pre - 
dicted by our models with the lReiprich fc Bohringeil J2002T) 
catalogue, again excluding those clusters with mass uncer- 
tainti es greater tha n 50 per cent, and using temperatures 
from iHornerl i200lf) . The shifted-entropy model provides a 
good description of the data over two orders of magnitude 
in mass, although it is not statistically the best fit. 

Unlike the M-L relation, the observations here show re- 
markably little scatter, and this scatter is consis tent with the 
publis hed measurement uncertainties. Recently. ISmith et alJ 
(2005) used strong gravitational lensing to measure the 
masses of ten X-ray luminous clusters and found that there 
is intrinsic scatter in the M-T relation, due mostly to merg- 
ing, non-equilibrium systems. Although this scatter is not 
apparent in Figure|S] it is possible that it is underrepresented 
here, as M200 is obtained from the X-ray data by assuming 
isothermality and i mposing a functio nal form for the surface 
brightness profile (iRasia et al]|2005l) . This parametrization 
could have the effect of homogenizing clusters with a range of 
cooling core sizes, geometries, and amount of substructure. 
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Figure 5. The relation between X-ray temperature and M2oo- 
The points are lo cal (2 < 0.2) data, with m asses from the HI- 
FLUGCS sample iReiprich &: Bohringe ill2002l) and t emperatures 
(uncorrected for cooling flows) from iHorner (2001) . 1— a error 
bars on the masses are shown; only galaxies with relative errors 
of < 50 per cent are plotted. The shaded region represents the 
prediction of the shifted-entropy model with K = 2i^ioo and a 
realistic 1<t (heavy shading) or 3<r (lighter shading) scatter in the 
halo concentration parameter. 



The same would also be true for the masses in the M-L re- 
lation, where we do see significant scatter (Figure |3J, and 
this could indicate that the intrinsic scatter in the luminosi- 
ties is larger than the intrinsic scatter in the temperatures. 
On the other hand, the fact that X-ray derived masses are 
directly proportional to the temperature introduces a cor- 
relation that could reduce the scatter in the M-T relation 
alone. It would be useful to have a larger sample of clus- 
ters with accurate lensing masses and X-ray observations to 
improve our understanding of the scatter in these relations. 

The model predictions in Figures |S] and HJ show that 
the predicted temperature is relatively insensitive to both 
the halo structure (i.e. concentration) and the entropy floor, 
for K > 0.5ifioo- Recall that the distribution of halo con- 
centrations partly arises from substructure, triaxiality and 
departures from equilibrium in simulated clusters and there- 
fore our predicted scatter approximately includes these ef- 
fects. However, in our model these concentrations are ap- 
plied to spherical, smooth halos and the scatter in pre- 
dicted temperatures is therefore not as large as in numerical 
and analytic models that do not make these assumptions 
l|Rowlev et alJl2004l : lOstriker et al.ll2005l) . The insensitivity 
of our predicted temperatures to the value of the entropy 
floor is due to the fact that the higher central tempera- 
ture associated with larger K a is offset by the flattening 
of the central density profile, which means the luminosity- 
weighted temperature is dominated by the temperature at 
larger radius. Thus, even a model with K a = 0. 5i\ioo pre- 
dicts temperatures that are just within the scatter of the 
observations. For the same reason, the K coo \ model predic- 
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Figure 6. As Figure but for the mass— temperature relation. 



tion does not differ significantly from those shown here, so 
we have omitted it from the figure for the sake of clarity. 

In summary, we have found that a variation of nearly 
an order of magnitude in K is required to explain the large 
scatter in the M — L relation. Variations in dark matter 
halo shape (concentration) alone are insufficient. Encour- 
agingly, the wide range of entropy levels required does not 
conflict with the small observed scatter in the M — T rela- 
tion, because of the temperature's inse nsitivity to the core 



properties. Although both theory (e.g. |Rrjwle^^^|j200 
lOstriker et alJ 120051) and observations Tsmffhet*^iTl20(T 



suggest that there may be additional scatter in the M-T 
relation that is not apparent under the simplifying assump- 
tions of both the models and data presented here, it is still 
less than the scatter in the M-L relation, which is dominated 
by variations in the entropy distribution of the gas. 



4 SCATTER AND THE LUMINOSITY AND 
TEMPERATURE FUNCTIONS 

If the mass spectrum of dark matter halos is known precisely, 
then the observed shape and normalization of the tempera- 
ture and luminosity functions provides an independent test 
of the theoretical models, that does not depend on an obser- 
vational determination of cluster mass. Scatter in the mean 
relations plays an important role here, and can influence the 
shape of these functions. 

We construct the theoretical luminosity and temper- 
ature functio n s usin g the dark matter mass function of 
Evj-arcL^^dJ tajCjd), based on the fitting formalism of 
Jenkins et alJ l|2001il , which provides a universal description 
of the mass fu nction to within ab out 10 per cent. The ad- 
vantage of the lEvrard et all l|2002f) mass function is that it 
is expressed in terms of M200, the same mass that we use 
to compare with the observed mass-temperature and mass- 
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Figure 7. The sha ded region is the ob served temperature func- 
tion at z = from llkebe et all |2002|) . bu t with t e mpera tures, 
uncorrected for cooling flo ws, ta ken from iHorneil feOOlD . The 
data points are from lHenrvl l2OO0t) . In the bottom panel we show 
the shifted-entropy model with three different lev els of he ating, 
assuming its = 0.85 as de termined from WMAP ISpergel et alJ 
1200.4 iTegmark et alj l20041. Models are convolved with a 10 per 
cent uncertainty in temperature which flattens the temperature 
function. In the top panel, we show the same models, normalized 
at 4-5 keV with different values of ag, as indicated. 



luminosity relations . The form of the temperature and lu- 
minosity functions are then completely determined by the 
correlation between virial mass and the X-ray observable 
(BBLP 3 ). 



4.1 The temperature function 

The observed z = 0.15 t e mper ature functions from 
llkebe et all <2002l) and iHenrvl (|200CT) are shown in Fig. U 
For co nsistency, the temperatures are taken from iHornerl 
(2001), and therefore uncorrected for any cooling flow com- 
ponent, though this makes little difference in practise. To 
compare with these data, we show predicted temperature 
functions from the shifted-entropy models with a range of 
normalizations that approximately accounts for the scatter 
in the M-L relation. We have smoothed the models with 
10 per cent Gaussian random noise on the temperatures, 
to mimic the scatter in the mean relation (which is consis- 
tent with being due to observational uncertainties). First, in 
the bottom panel, we show three models with different K a 

2 We note that M500 is better determined observationally, and it 
would be useful to have a mass function from numerical simula- 
tions filtered on this scale. 

3 However, both the observed and model luminosity functions 
shown in Figure 9 of BBLP are incorrect due to errors in the 
bolometric correction and cosmology conversion. 
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but the same normalization as — 0.85, as measured from the 
WMAP data (and adjusted for the slightly non-concordance 
values of cosmoloK i cal parameters that we have adopted 
ISpereel et ail 120031: Freemark et all 12004) . At the hot end 
of the temperature function, the data is best matched by 
the models with the highest entropy floors, while at the op- 
posite extreme the low-entropy models fare better. However, 
we caution that at low temperatures (< 1 keV) there may 
be completeness issues that could artificially flatten the tem- 
perature function and, therefore, yield values of K a that are 
systematically lower than that of the average system (e.g. 
lOsmond fc PonmadEloll) . 

In the top panel of Fig. |7| we again show the predicted 
temperature functions for a range of K a values, but with 
as adjusted to give the same number density of clusters at 
T — 4-5 keV, where observ ational data from different stud- 
ies are in best agreement dlkebe e t al. 2002). The best-fit 
value of as is ~ 3.5 per cent higher if the normalization to 
the data is made over the range kT — 6-7 keV. The range 
of K (which are all reasonably consistent with the mass- 
temperature relation) corresponds to a range of best-fit val- 
ues of as ranging from 0.76 to 0.9. This is com parable to the 
obser vational uncertainty on this parameter dSpergel et alJ 
2003), and therefore we cannot use the temperature function 
alone to provide a sensitive test of the size of the entropy 
floor. This is again simply because temperature is relatively 
insensitive to the entropy of the central gas. 

4.2 The Luminosity Function 

To compute the luminosity function, the intrinsic scatter in 
the M-L relation must be taken into account. Because of the 
steepness of the mass function, even a small distribution of 
halo masses corresponding to a given luminosity can have 
an important effect on the number density of clusters at 
that luminosity. Unfortunately, the value of K in our model 
does not have a unique physical motivation, and thus we do 
not have a prediction for the scatter as a function of mass. 
However, we can see from the data in Figure |H] that the 
observations are approximately covered by models with a 
range of entropy floors 0.5 < K /Kioo < 7, so in Figure|H|we 
show the prediction of the luminosity function for these two 
extremes. The true luminosity function should lie between 
these limits, with a shape that depends on the distribution 
of entropy levels at each luminosity. 

The obse rved luminosity f unctions at z = 0.15 from 
the WAR PS llJones et al 1 120011 Jones et al., in prep.) and 
REFLEX jBohrineer et alJl200ll) surveys are shown in Fig- 
ure |Sj The solid lines show our model, for K a = 0.5-Kioo 
(upper line) and K = 7K\oo (lower line). This range brack- 
ets the observational data, although the data do lie nearer 
the model with high entropy. This is especially true for the 
low-luminosity clusters, with L % 10 44 ergs s _1 . This may 
indicate that lower-mass clusters have higher central en- 
tropies, on average; this may also be evident from Figure [3] 
although there are few clusters with accurate mass measure- 
ments at these low luminosities. 

We also show, as the dotted lines, the effect of vary- 
ing as between 0.76 and 0.9 (assuming K — 2_Kioo). This 
has only a small effect on the bright end of the luminosity 
function, and no effect on the faint end. The dashed line 
shows the effect of removing the limiting surface brightness 
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Figure 8. The lines in each panel show different theoretical mod- 
els for the luminosity function. The solid lines show models with 
entropy floors of K = O.bKiQQ (upper line) and K a = 7-K"ioo 
(lower line). This approximately brackets the range of entropies 
required to explain the scatter in the M-L relation (Figure^- The 
dotted lines show the effect of varying erg as indicated, keeping 
K = 2K\qq fixed; lower values of as reduce the number of lumi- 
nous clusters by a small amount but have no effect on the faint 
end of the luminosity function. Finally, the dashed line shows the 
K = 2A"ioo model but omitting the surface brightness thresh- 
old, which increases the prediction at low luminosities. Left panel: 
The open squares with error bars are the observed local luminos - 
ity function from the REFLEX survey iBohringer et alJl20o3) . 
Right panel: Simil ar, but where the data are from the WARPS 
Jjones et al.ll200lL Jones et al., in prep.), in a different energy 
band. 

thres hold of 1 x 10~ 15 ergs s -1 cm -2 arcmin -2 iScharf et alJ 
1997) used in the other models; these surface brightness cor- 
rections are relevant only for the least luminous clusters in 
the sample. 

Thus we have shown that consistency between the X- 
ray scaling relations (M-L and M-T) and the luminosity 
function can be achieved in these models; however, a better 
understanding of the entropy-floor distribution as a func- 
tion of mass is required to make a firm prediction of the 
luminosity function shape. 



5 EVOLUTION 

We now turn to the redshift evolution of the X-ray scaling re- 
lations and the temperature and luminosity functions. These 
predictions can provide another interesting test of the differ- 
ence between the fixed-i^o, preheating models and models 
where K a = K coo \ . 

In Fig. |5| we show the predicted evolution in the mass- 
temperature and mass-luminosity relations for the two mod- 
els. The data are the same z ~ 0.15 data shown in Fig- 
ures and and the model predictions are shown at 
z = 0,0.15,0.4 and 0.7. For massive clusters, T £ 4 keV, 
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Figure 9. The observed M— L and M— T relations shown in Figs, 131 
and |S| are reproduced as the solid circles. In the bottom panels 
we show the default model predictions, where the entropy floor is 
independent of redshift, at z = (solid line) and z = 0.15, 0.4, 0.7 
(dotted lines). The curves in the top panels are the models where 
the entropy floor K coo i is related to the cooling time of the gas 
and thus evolves with redshift. 



the predicted evolution in both the M-T and M-L correla- 
tions is mostly in the normalization, with little change in 
the slope. 

The amount of evolution in the M-L relation is small 
relative to the observed scatter at z = 0. On the other hand 
predicted evolution in the M-T relation is more noticeable; 
clusters at a given temperature are predicted to be 40-60 
per cent less massive at z ~ 0.7. The sense and magnitude 
of the evolution are comparable to recent XMM-Newton 
and C handra data llKotov fc Vikhlininl2005l : lMaughan et alJ 
l2006h . Both the K a = 2K W0 and the K a = K coo i models 
models predict a similar amount of evolution for the M — T 
relation, so this is not a useful way to discriminate between 
them. 

Interestingly, although the evolution in the M-T and 
M-L relations appear similar for both models, the predicted 
evolution of the L-T scaling law is in opposite directions, 
as shown in Figure 1101 At high temperatures, the amount 
of evolution in the predicted relation is very small. For the 
fixed-floor model, the evolution is negligible, while the K coo \ 
model predicts that high redshift clusters will be about 
30 per cent more luminous at fixed temperature, due to 
the fact that the entropy floor is lower at higher redshift. 
Observations of distant clusters seem to indicate a much 
stronger evolution, with clusters at z ~ 0.7 being up to 2.5 
times brighter than local clusters of the same temperature 
llVikh lmm et al.ll2002l:lLumb et al.ll200 4: Kot ov fc Vikh linin 
120051 : iMauehan et alJl2006T) . In principle, this observation 
has the potential to rule out both models shown here. How- 
ever, we note that the amount of evol ution observed is sen si- 
tive to the details of the analysis (e.g. lEttori et alJl2004allbl) . 
Since any evolution in the mean scaling relation is much less 
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Figure 10. The observed temperature-luminosity relation at z = 
is shown as the solid circles. The thin and thick lines represent 
the K coo i and K a = 2i\ioo models, respectively. Predictions are 
shown for z = (solid line) and z = 0.15,0.4,0.7 (dotted lines). 

than the factor ~ 4 scatter in luminosity at fixed temper- 
ature for local clusters, it is probably premature to claim 
these observations rule out either model until the scatter 
and the selection biases that result from it (i.e. Malmquist- 
like bias) are robustly integrated into the model predictions. 
The most leverage will come from low temperature systems 
at high redshift; the K a = 2iTioo model predicts very strong 
(negative) evolution in the luminosities of these groups, as 
the fixed entropy floor becomes very large relative to the 
characteristic entropy. 

The evolution of the temperature and luminosity func- 
tions requires a knowledge of the mass function at z > 0. 
Since this has not yet been precisely measur ed from simu- 
lation s, we take the local mass function from lEvrard et ail 
(2002) and evolve it to higher redshift using Press-Schechter 
theory. In Figure [TTI we show the temperature function for 
both models (normalized at kT ~4-5 keV, with erg ~ 0.8), 
at z = 0, 0.4 and z — 0.7, compared with z = 0.4 data from 
iHenrvl (2000) 4 . Over the range of observed temperatures, 
kT > 3 keV, the evolution in the two models is similar, and 
in good agreement with the observations. The models begin 
to diverge at lower temperatures, where the fixed-entropy 
model predicts a little less evolution. 

The predicted evolution of the luminosity function is 
shown in Figure H2l Unfortunately, our model does not pre- 
dict the scatter in the M-L relation nor its evolution, so 
for illustration we have just shown the K a — 27fioo model, 
which provides a reasonable match to the bright end of the 
local luminosity function. The models are again compared 
with data at sim ilar redshifts, from the WARPS survey 
Jjones et al.ll200ll Jones et al, in prep.). The WARPS data 

4 An updated version of the observed high-r edshift t empe rature 
function, based on more data, is presented in IHenrvl (2004). The 
results are consistent with those shown here. 
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Figure 11. The predicted temperature functions at z = 0,0.4 
and z = 0.7, for the fixed— entropy models (thick lines) a nd the 
A" nnn i models (thin lines). The data are at z ~ 0.4, from lifenrvl 
(2000). 

at z>0.6 have been updated with accurate luminosities mea- 
sured from XMM-Newton and Chandra observations, and 
are corrected to the rest-frame energy range 0.5-2.0 keV. 
In contrast with the temperature function, the observed lu- 
minosity function evolution is modest, with a factor ~ 3 
decrease in the number of massive clusters between z — 
and z = 0.7. 

The ifcool model predicts very little evolution in the 
luminosity function, and thus overpredicts the number of 
high redshift clusters. The K a = "2K\qo model appears to 
be in much better agreement with the data, especially for 
the brightest clusters. As with the local luminosity function, 
the disagreement at the faint end may indicate that lower 
luminosity clusters have larger central entropies. However, 
as we have already seen (Figure 0, the shape of the lumi- 
nosity function is very sensitive to the amount of scatter 
in the M-L relation, and we have no theoretical or empiri- 
cal knowledge about how this scatter evolves. We also note 
that the high-luminosity end of the luminosity function is 
still poorly determined, with a variation in observed num- 
ber abundance at fixed luminosity measured f rom different 
surve ys at z > 0.3 being about a factor ~ 2 (Mulli s~et alJ 
2004). 



6 DISCUSSION AND CONCLUSIONS 

In this paper, we have revisited the observed X-ray scal- 
ing relationships between mass and temperature (M-T) and 
mass and luminosity (M-L), simultaneously with the tem- 
perature and luminosity functions. We have attempted to 
compare the observable quantities as directly as possible (i.e. 
without bolometric or cooling-flow corrections) with a sim- 
ple suite of models in which a fiducial gas entropy distribu- 
tion is shifted by a value 0.5 < K a /Kioa < 7. In particular, 
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Figure 12. The luminosity functions at z = 0,0.4 and z = 0.7, 
for the fixed— entropy (2i^ioo) models (thick lines) and the A" coo i 
models (thin lines). Data at z = 0.4 and 0.7 are shown, from 
the WARPS survey, updated with XMM-Newton and Chandra 
luminosities for clusters at z > 0.6. 

we focus on the scatter in the observed scaling relations, 
and how this compares with the scatter expected due to 
a) a range of halo structures (concentrations); b) the time 
available for cooling-only processes or c) heating/cooling ef- 
ficiency. Since cluster temperatures are relatively insensitive 
to variations in the entropy distribution (and hence the scat- 
ter in the M-T relation is small), we gain the most by focus- 
ing on the M-L relation and the luminosity function. Our 
main findings are as follows: 

(i) The variations in dark matter halo concentration ex- 
pected from simulations are not large enough to account 
for the scatter in the observed M-L relation of clusters and 
groups. 

(ii) Simple models of the intracluster medium in which 
the core entropy is modified to have a minimum value require 
the value of this floor to be between about 0. 5-Kioo and 
7-Kioo to match the slope, normalization and scatter in the 
observed M-T and M-L scaling relations. The constraint 
comes mostly from the M-L relation, as the temperature is 
insensitive to the value of K a . 

(iii) The shape of the luminosity function is sensitive to 
the scatter in the M-L relation. The observations lie between 
the models with K a = 0.5ifioo and K = 7A"ioo, but closer 
to the higher-entropy model. The scatter in entropy levels 
as a function of halo mass must be accounted for if the 
parameters as or K a are to be accurately deduced from the 
luminosity function alone. 

(iv) The model temperatures are in good agreement with 
the obs erved tempera t ure fu nction, assuming the mass func- 
tion of lEvrard et al.l B2002I) . However, the insensitivity of 
temperature to K , and the uncertainty on the normaliza- 
tion parameter as, means this does not put strong con- 
straints on the value (or range of values) of K . 

(v) The amount of gas that can cool in a Hubble time sets 
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a maximum value on the minimum entropy of intracluster 
gas, K - Lower entropy floors can be achieved by allowing 
some gas to cool below the threshold, or reducing the time 
available for cooling. However, although this model provides 
a reasonable match to the median mass-luminosity relation, 
it cannot account for the many clusters with luminosities 
below this relation. 

(vi) We also present predictions for the evolution of the 
scaling relations and temperature/luminosity functions, out 
to z = 0.7. A comparison with observations suggests that the 
model with K a = 2.K100 is valid out to z ~ 0.7, but this de- 
pends on how the scatter in the M-L relation evolves, which 
is currently unknown. On group scales, where observations 
unfortunately will be most difficult, these predictions clearly 
distinguish between an entropy floor K that is independent 
of mass and redshift, and one that is tied to the cooling time 
of a halo. 

Thus, the scatter in the observations can be understood 
if the minimum entropy of the gas in clusters has a median 
value of K = 200 keV cm 2 but varies by a factor ~ 3 
between halos of similar mass. Although we lack a quanti- 
tative theoretical prediction for the origin of this scatter, it 
seems reasonable to associate it with a similar range in the 
efficiency of heating and/or cooling. 

We acknowledge that the models presented here pro- 
vide an incomplete description of the intracluster medium. 
Apart from the obvious point that the pure heating mod- 
els neglect the cooling processes that must take place to 
form the galaxies and AGN responsible for the heating in 
the first place, these models also lead to isentropic core gas 
distributions that are in conflict with a t least some high res- 
olution observations o f clusters (e.g., iPonman et al.ll2003t 
iPratt fc Arnaudll2005t) . Therefore, cooling must be incorpo- 
rated for a full description of the data. Although we have 
made a very crude step in this direction by tying the en- 
tropy floor level to the cooling time available, this model is 
greatly oversimplified, as it does not allow the gas to flow 
to the centre as it cools as in more realistic models (e.g. 
VBBB). Furthermore, for low mass clusters and groups this 
model predicts too much condensed g ELS, EIS most of the in- 
tracluster medium in a self-similar model has an entropy 
below the cooling threshold. Therefore, some combination 
of heating and cooling is expected to be required. 

Such a model has recently been developed by 
iMcCarthv et alJ (l2004fl : in this model, preheated gas is al- 
lowed to cool in a realistic way, following the hydrodynamic 
evolution of the gas as it flows to the centre. With an ap- 
propriate range of preheating levels and cooling times, this 
model can reproduce the slope, normalization and the scat- 
ter in the M-L relation. Our results show that including 
scatter in the halo potentials (not considered by McCarthy 
et al.) would only have a secondary impact on their con- 
straints. We have also explored the effect of scatter on the 
temperature and luminosity functions. Since the M-T rela- 
tion is relatively insensitive to the gas distribution (as long 
as the very lowest entropy gas is heated or removed), we can 
expect that the predictions for the temperature function us- 
ing the model of McCarthy et al. would be very similar to 
the models presented here. However, the luminosity func- 
tion, and the evolution in the scaling relations, will be sen- 
sitive to the mass and redshift dependence of the distribu- 



tions of preheating levels and cooling times. The next step, 
therefore, is to physically link these distributions through 
a feedback model, and to compare the resulting luminosity 
function with the observations. 

Finally, we note that while analyses of the M-L and 
M-T relations and the luminosity and temperature func- 
tions can teach us much, they are only probing the inte- 
grated properties of the intracluster medium. Some of the 
most powerful constraints on the role of non-gravitational 
physics in mediating the properties of the ICM are likely 
to come from detailed comparisons of theoretical models 
to the actual spatially-resolved profiles (e.g., entropy and 
temperature) of clusters. The difficulty is in obtaining such 
profiles for large, representative samples of clusters and (es- 
pecially) groups. However, as more and more groups and 
clusters observed with Chandra and XMM-Newton become 
publicly availa ble, progress is beginning to be made on 
this front (e.g. IPratt fc Arnau dll2005l : iPiffaretti et"al"ll2005l ; 
ICvpriano et alJl200Sft . Such comparisons will make excel- 
lent complementary probes to studies such as the present 
one that make use of more easily obtained and robust inte- 
grated cluster properties. 
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